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Abstract 



Quarkonium production in hadron collisions at low transverse momentum can be used for probing 
Transverse Momentum Dependent(TMD) gluon distributions. For this purpose one needs to establish 
the TMD factorization for the process. We examine the factorization at one-loop level for the produc- 
tion of r\ c or r]b. The perturbative coefficient in the factorization is determined at one-loop accuracy. 
Comparing the factorization derived at tree-level and that beyond the tree-level, a soft factor is in 
^ c"| general needed to completely cancel soft divergences. We have also discussed possible complications 

| of TMD factorization of p-wave quarkonium production. 

>: 

Quarkonium production in hadron collision can be used to explore the gluon content of hadrons, 
because the quarkonium is dominantly produced through gluon-gluon fusions. For the produced quarko- 
nium with large transverse momentum, one can apply QCD collinear factorizations for long distance 
effects of the initial hadrons. In this case, one can extract the standard gluon distributions (see e.g., [1]). 
If the quakonium is produced with small transverse momentum q±, it can be thought that the small q± 
is generated at least partly from the transverse motion of gluons inside the initial hadrons. In this case 
one can apply Transverse Momentum Dependent (TMD) factorization for initial hadrons. Therefore, the 
production with small q± allows to access TMD gluon distributions. 

Factorizations with TMD quark distributions and fragmentation functions have been studied inten- 
sively beyond tree-level in different processes in[2l [3J, HI [5]. In comparison, the factorization with TMD 
gluon distributions beyond tree-level has only been studied for Higgs production in hadron collision in 
[6]. Recently, the TMD factorization of quarkonium production has been derived at tree- level in [7], and 
based on it numerical predictions have been obtained. For theoretical consistency and precision it is 
important to examine the TMD factorization beyond tree-level. From early studies in [2j (3] [H [6] it is 
known that a soft factor needs to be implemented into the factorization. In this work we examine TMD 
factorization of r/ c or % production at one-loop level. 

A quarkonium is dominantly a bound state of a heavy quark Q and its antiquark Q. Because of the 
heavy mass the QQ pair is of a nonrelativistic system. To separate the nonperturbative effects related to 
the quarkonium in its production one can employ nonrelativistic QCD(NRQCD) factorization [8] by an 
expansion of the small velocity of Q relative to Q . The inclusive production of a quarkonium at moderate- 
or large q± has been studied intensively both in theory and in experiment. In last five years, important 
progresses were made in the study of the next-to- leading order QCD correction for J/ip production in 
hadron collisions [9 J and power corrections [10J. The activities in this field can be seen in It should 
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be noted that in experiment it is also possible to study the inclusive production at low qj_. E.g., a 
J/i/j produced at LHCb can be measured with q± smaller than !GeV[T2]. Therefore, with theoretically 
established TMD factorization one can extract from experimental results TMD gluon distributions. 

We will use the light-cone coordinate system, in which a vector is expressed as = (a + ,a~ , aj_) = 
((a + a 3 )/\/2, (a - a 3 )/\/2, a 1 , a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 . We introduce two light cone vectors = 
(0, 1, 0, 0) and = (1, 0, 0, 0) and the transverse metric = g^ u — n^l 1 ' — n^l 11 . We consider the process: 

h A (P A )+h B (P B )^ VQ (q) + X, (1) 

in the kinematical region Q 2 = q 2 3> q\ with Q as the mass of rjq, where tjq stands for r/ c or 775. The 
momenta of the initial hadrons and of the quarkonium are given by 

P^ (P+ 0,0,0), Pgw (0,^,0,0), q^ = (xP+,yP B ,q ± ), (2) 

where we have neglected masses of hadrons, i.e., ~ and Pg ~ 0. For each hadron in the initial 
state one can define its TMD gluon distribution. We introduce the gauge link along the direction = 
(u + ,u~, 0,0): 

£ u (z, — 00) = Pexp ^— ig s J dXu ■ G(Xu + z)\ , (3) 

where the gluon field is in the adjoint representation. At leading twist one can define two TMD gluon 
distributions through the gluon density matrix [13]: 



xP+ J (2vr) 

= -29±fg/A(x,k±,£,v) + + \a± k \) h g / A (.x,k ± ,C 2 ,n) (4) 



with = (0,^ _ ,^j_). x is the momentum fraction carried by the gluon inside h A . The gluon has also a 
nonzero transverse momentum k±. The definition is given in non-singular gauges. It is gauge invariant. 
In singular gauges, one needs to add gauge links along transverse direction at £~ = — oo|14|. Because of 
the gauge links, the TMD gluon distributions also depend on the vector u through the variable: 

In the definition the limit u + <C u~ is taken in the sense that one neglects all contributions suppressed 
by negative powers of Q 2 . 

From the definition in Eq.@ there are two TMD gluon distributions. The distribution f g / A cor- 
responds to the standard gluon distribution in collinear factorization. The distribution h g / A describes 
gluons with linear polarization inside h A . The relevant phenomenology of h g / A has been only recently 
studied [151 [16]. Through the process studied here one can also obtain information about this distribution 
[7]. For \ib one can also define two TMD gluon distributions f g m and h g / B similar to those in Eq.@, 
in which the gauge links are along the direction = (v + ,v~ ,0,0) instead of u and the limit v + 3> v~ 
is taken. Therefore the two distributions f g m and h g /B depend on the parameter ( v which is defined by 
replacing in C, u P A with Pg and u with v in Eq.([5]). 

To study the TMD factorization of the process in Eq. ([1]) , we need to study 

g(p)+g(p)^r] Q (q)+X, (6) 
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with p^ = (Pti 0,0,0) and p^ = (0,P^",0, 0). The reason is the following: If the factorization holds, it 
holds for arbitrary hadrons in the initial state. Especially, it also holds if the initial states are of partons. 
Correspondingly, we need to study the TMD gluon distributions of the corresponding gluon. That is, we 
need to study f g /A and h g / A of the gluon g{p) instead of \ia and f g /g and h g /B of the gluon g(p) instead 
of Kb • At tree-level, one easily finds: 



k±, (i, m) = f^/ B {x, k ± ,C,v) = 8(1 - x)5 2 (k ± ), 



r(0) 



(7) 



while h g /A and h g m are zero. They become nonzero at order of a s . 

In this work we will work at the leading order of the small velocity expansion in NRQCD. At this 
order the production of tjq can be thought as a two-step process. In the first step a QQ pair is produced in 
which the heavy quark Q and its antiquark Q carry the same momentum q/2. The pair is in color-singlet 
and spin-sing let 1 S - Then the pair is transmitted into tjq with the mass Q = 2toq. The transition is 
described by an NRQCD matrix element. At tree-level, the process in Eq.([6]) is with X as nothing. It is 
straightforward to obtain the differential cross-section: 



da 



7T 



dxdyd 2 qx_ 



cr -^S(xys - Q 2 )5(l - x)S(l - y)S 2 (q ± ), 



(AlTOis 



N C (N 2 - l)m Q 



(8) 



with s = 2p + p~ and mg being the pole mass of the heavy quark. V(0) is the wave function of t)q at the 
origin. In fact ^(O)] 2 should be expressed as a NRQCD matrix element. With results in Eq.([7|) one can 
derive the factorization: 



da 



dxdyd 2 q±_ 

n 



TTOq 

Q 2 
i + %) 



d 2 k a± d 2 k b± f g/A (x, k a ±)fg/ B (y, h±)6 2 (k al _ + k bl _ - q±)5(xys - Q 2 )U, 



(9) 



at tree-level. Beyond the tree-level one needs to introduce a soft factor. As we will see explicitly , all soft 
divergences will be factorized into the soft factor and TMD gluon distributions so that the perturbative 
coefficient % is free from soft divergence. We will then determine T~L at one-loop level. 
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+ h.c. 
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Figure 1: The one- loop corrections to the gluon TMD. The double lines represent the gauge link. The 
black bubble in Fig. la is for self-energy correction. 

To derive the factorization at one-loop we need to study the one-loop corrections to TMD gluon 
distributions and the differential cross-section. The one-loop correction to TMD gluon distribution has 
been studied in [6] where the collinear divergence has been regularized with a infinitely small off-shellness 
of the gluon. Here we regularize all divergences in d = 4 — e space-time. The correction can be divided 
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into the virtual- and real correction. The virtual correction is given by diagrams in Fig.l. We will use 
MS-scheme to subtract U.V. divergences. After the subtraction we have the virtual correction from Fig.l: 



fg/Ai X l ^J-J CliJ A 4 ) 




5tt z 
~Y2 



e s An/jj 



(10) 



where the poles in e s = 4 — d stand for collinear- or infrared divergences, i.e., soft divergences. fj, s is the 
scale associated with these poles. \x is the U.V. scale. The terms in the first line in Eq.(10) is the sum of 
the contributions from Fig. la and Fig.lc with their conjugated diagrams. The remaining terms are from 
Fig. lb and its conjugated diagram. 
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Figure 2: The real correction at one- loop to the gluon TMD. The double lines represent the gauge link. 
These diagrams are for real corrections. 



The corrections from Fig. 2 are real corrections. They can be found in [6] as 



a s N c 
ir 2 k 2 



x 



+ 



(1 



+ x(l — x) + 



X 1 r/i 

- + -5(1 
2 2 v 



1 



5(1 -x) 



x) In 



c 2 

k 2 



(11) 



where the terms in the first line are from Fig. 2a and Fig. 2d. The total one- loop correction is then the 
sum of the virtual- and real correction. At one-loop h g i^ becomes nonzero. It receives a contribution 
from Fig. 2a. We have: 

2a s N c 1 — x 2 

By replacing ( u with ( v we obtain f g i B and h g m from f g iA an d h g /A, respectively. 

Now we turn to one-loop corrections of the differential cross-section. The corrections can be divided 
into the virtual correction and the real correction. The virtual correction is the one-loop correction to 
the process g(p) + g(p) — > rjQ(q). We denote the total contribution from the virtual correction as 



do-(gg -)• t)q) 



dxdyd 2 q± 



2s(2vr) : 



Mxys - Q 2 )5(l - x)5(l - y)5 2 (q ± )a 1 . 



(13) 
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Figure 3: The class of diagrams where a gluon is emitted from the initial gluon g(p) and is attached to 
a possible place. There are 6 diagrams. 3 of them are given here. Another 3 diagrams are obtained by 
reversing the direction of the heavy quark line. 



The contributions to <j\ can be divided into 4 parts: 



<?\ = <?\A + <t\b + eric + &iD, 



(14) 



(T\A receives contributions from diagrams in which a virtual gluon is emitted by the initial gluon g{p). 
The diagrams for this part are given in Fig. 3. <j\b receives contributions from diagrams in which a virtual 
gluon is emitted by the initial gluon g(p). o\c denotes the contributions from diagrams in which a virtual 
gluon is exchanged between heavy quark line. a\o denotes the one-loop corrections of external gluon 
lines. This part will not contribute to H, because the contributions to <t\d are automatically subtracted 
into TMD gluon distributions. In the below we will only give and discuss the results of (t\a,ib,iC- 

In the above classification Fig. 3a can contribute both to ct\a and o\b- We put the half of the 
contribution Fig. 3a into o\a and another half into o\b- With symmetry arguments one easily finds 
&1A = a lB- We have then: 



1 

&XA = &1B = 2°"1 





+ <Tl 




3a 


3b 


3c 



(15) 



By summing contributions from each diagram we obtain the following results for the virtual corrections: 

,2 
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(TO 



OsNc 

12vr 



4 2 / e^ATiul 
6--6-(l + ln / s 



3 In 2 - 6 In ■ ' 



Q- 



Q 2 



ri 2 11 
+91n^ -61n2 + 6 + — 7T 
Q 2 4 

2 



2vr 



iV c ln^ + C7 F ^-2 + 41n2j +-L(^21n2-^ 2 



(16) 



In these results the U.V. poles are subtracted in the MS-scheme. The on-shell scheme for the renormal- 
ization of heavy quark propagators is used so that rriQ = Q/2 is the pole mass of heavy quark. In g\a 
the pole terms of e s are for soft divergences coming only from Fig. 3a. The contributions from Fig. 3b and 
Fig. 3c also contain collinear divergences and infrared divergences. The infrared divergences are canceled 
in the sum of the two diagrams, because the QQ is in color singlet. The collinear divergences are also 
canceled. In calculating the diagrams for o\c one will meet Coulomb singularity. This singularity is 
factorized into NRQCD matrix element. Hence we have finite o\c- 
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The real correction is from the tree-level process 

g(p) + g(p) ->■ r) Q (q) + g(k). (17) 

For the color-single QQ-p&ir there are 12 diagrams for the amplitude. Since we are interested in the low q± 
region, we expand the differential cross section in q±/Q and only take the leading order in the expansion. 
At the leading order we have only those diagrams given in Fig. 4 for the differential cross-section. The 
result for the process in Eq. (jl7j) in the limit of q±_ — > is: 




Figure 4: The diagrams for the cross-section of g + g — > g + t]q . In these diagrams, the gluon in the 
intermediate state is emitted or absorbted by gluons. The black dots denote the projection of the QQ 
pair into the color singlet 1 So state. By resversing the quark lines one can obtain other 3 diagrams from 
each diagram. 

The factorized result in Eq.Q is derived at tree- level. If we extent the factorization beyond tree- 
level, with the one- loop results in the above we will find the following: 1). The soft divergences are 
not factorized, i.e., % will contain some infrared divergences represented by poles in e s . 2). The real 
correction of the differential cross-section is not totally generated by TMD gluon distributions, in other 
word, T~L will receives correction from the real correction. It results in that % depends on q±. All of these 
have the common reason. In the one-loop corrections to the differential cross-section there are exchange of 
a soft gluon between the two initial gluons g(p) and g(p). In the virtual correction the exchange results in 
infrared divergences, and in the real correction it results in contributions proportional to 5(1 — x)S(l — y). 
The effects of the soft gluon exchange are not exactly generated by the corresponding soft gluon exchange 
in TMD gluon distributions. The effects of soft gluon exchange are of long-distance. Therefore, one needs 
to introduce a soft factor in the factorization to completely factorize these effects from % determined 
with Eq.@. 

The effects of soft gluon exchange between a gluon moving in the +- direction and a gluon moving 
in the — direction can be described by the expectation value of a product with four gauge links. We 
introduce as in [6]: 

S(b±,fA, p) = ^2^7<0|Tr [4(6±,-oo)£ M (6 ± ,-oo)4(o,-oo)£„(0,-oo)] |0>. (19) 

The gauge links are past-pointing. It reflects the fact that the two gluons g(p) and g(p) are in the 
initial state. The dependence on the directions of gauge links is only through the parameter p 2 = 
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(2u • v) 2 /(u 2 v 2 ) ~ u~v + /(u + v~). The limit u~ <C ii + and v + S> w~ is taken similarly to that in TMD 
gluon distributions. The gauge links or the gauge field is in the adjoint representation. At leading order 
one has 

S^(b ± ,p,,p) = l. (20) 



+ h.c. 




Figure 5: One-loop corrections for the soft factor. The first three diagrams plus their complex conjugated 
are virtual corrections. The last four diagrams are real corrections. A cut line is implied. 

At one-loop there are corrections from Fig. 5. One can divide the corrections into a virtual- and a real 
part. The diagrams in the first row are of the virtual part. Those in the second row are of the real part. 
The virtual correction read: 



u s N c 
2vr 



e s 4vr/i2 



(2 -In/) 



(21) 



where the U.V. pole is subtracted. The pole in e s represents the I.R. divergence with the scale p s . The 
real part is: 



sW(b ± ^,p) = -^(2-ln P 2 ) 



2vr 2 



d 2 k,- 



k\ 



(22) 



The total one-loop contribution is the sum of the virtual- and real contribution. 
We now define our soft factor which will enter the TMD factorization as: 



S(£±,p,p) 



d 2 b ± 
(2vr) 2 



a s N c 
2ir 



(2 -hip 2 ) 



2 | h eV 
e s Airp 2 



6 2 ( 



+ 0(a 2 ) 



(23) 



With the introduced soft factor we propose the TMD factorization as: 



da 



dxdyd 2 q± 



TTOQ 

Q 2 



^ J d 2 k a± d 2 k b± d 2 i ± 5 2 (k a± + k b± + £ ± - q ± )5(xys - Q 2 ) 

■fg/A(x, k a _l_,(u, P)fg/B(y, h±,Cv,^)S(£±,p, p)n(Q, p, Cu,Cv)- 



(24) 
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Prom one-loop results of the differential cross-section, TMD gluon distributions and the soft factor we 
derive: 



H{Q, p, Cu, Cv 



1 + 



a,N c 



47T 



r 2 c 2 

l n 2 W. + l n 2 iw. 

Q 2 + Q 2 



lnp 2 l + 21n 



Q 2 



+ 



iV 2 



1 



-7T 



+ 0{a 2 s ). 



+ 21n^ + ^ 



(25) 



It is clear that T~L is free from any soft-divergence and does not depend on q±. With the factorization 
the small transverse momentum q± is generated by the transverse motion of gluons in the initial hadrons 
and by soft gluon radiation. Eq. (124p and Eq. (j25j) are of our main results. It should be noted that 
the factorization holds for arbitrary large ( u and £ v . For practical applications one may take a frame 
to simplify the results in Eq. (|24p and Eq. (|25p . One can take Cu = Cv = PQ 2 so that the ™ D g luon 
distributions in Eq. (|24p depends on p and Q 2 and the perturbative coefficient becomes a function of Q, 
p and p: 



n(Q,p, P ) = i + 



a s N c 
2vr 



ln 2 p 



In p 



1 + 21 ^ 



/i 



7 



1 



+ In + -tT + — r 



Q 



N 2 

c 



1 , 



(26) 




Figure 6: The diagrams for the cross-section of g + g — > g + g + t]q give the contributions factorized with 
the gluon TMD g^. By resversing the quark lines one can obtain other 3 diagrams from each diagram. 

At the orders we have consider we do not find the contributions which are factorized with h g /^ or 
hg/g defined in Eq.@. But, there is a contribution involving these distributions of linearly polarized 
gluons in the TMD factorization. This contribution can be found at higher order of a s from diagrams 
given in Fig. 6. It is straightforward to calculate these diagrams in the limit q± —> 0. We find that the 
contribution takes the factorized form: 



da 



dxdyd 2 q± 



Fig.6 



^-5(xys - Q 2 ) I <f 2 l,„ . d-h,,() 2 [l + k h _ - (/_]_) 



fg/Ai x i h 



\{$a±. ■ h±f - ^a± fc 6± )h g / A (x,k al _)h g/B (y,k b j_) 



2a 



f g /B(y,h±] 



2a 

(27) 



Our result in the last line has also been derived in [7] with a different method. The perturbative coefficient 
of the contribution in the last line is at order of a°. This contribution should be added to Eq. (l24p . In 
principle one can determine the perturbative coefficient of the contribution beyond the leading order of 
a s following the same way as has been done for Eq. (|24|25p . But this will be very tedious because one 
needs to calculated the partonic process g + g — > t}q + X at 3-loop level. We leave this for future study. 
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The studied TMD factorization can be used for the region with qj_ ~ h-QCD for extracting TMD 
gluon distributions. But its usage is not limited to this kinematic region, because the factorization holds 
in general in the region gx/Q « 1. In the region Q > 9i > ^qcd both TMD factorization and 
collinear factorization hold. In the collinear factorization the perturbative coefficient functions in this 
region contain large log of qj_/Q- The results from TMD factorization can be used to resume these large 
log's. This leads to the well-known CSS resumation[3]. Based on our result here one can also derive the 
resummation in the case of quarkonium production, similarly to that derived in [6]. We therefore do not 
discuss the details about the resummation here. We note that such a resummation has been studied very 
recently in [T7j. An early work about the resummation can be found in |18j . where the formation of a 
quarkonium from a QQ pair is described with a color evaporation model instead of NRQCD factorization. 

At the orders we have considered, the production of a p-wave quarkonium is possible. But, the TMD 
factorization in this case can be complicated. According to the NRQCD factorization in [8], one needs 
to consider not only the contribution from the production of a color singlet p-wave QQ pair, but also 
the contribution of a color-octet s-wave QQ pair. The formation of a p-wave quarkonium from the color- 
singlet- and the color-octet QQ pair is at the same order in the small velocity expansion. In the case we 
studied here, we only need to consider the contribution from production of a color-singlet s-wave QQ pair. 
At the leading power the pair decouples with soft gluons. However in the case of p-wave quarkonia, the 
color-singlet p-wave- and color-octet s-wave QQ pair can emit soft gluons at leading power. To completely 
separate the effects of soft gluons, one may need a different soft factor than that introduced here. This 
is also the reason why we write our TMD factorization in Eq. (|24p explicitly with the unsubtracted TMD 
gluon distributions and the soft factor. It is noted that one may re-define TMD gluon distributions to 
factorize the differential cross section only with the re-defined TMD gluon distributions [19J. Another 
complication with p-wave quarkonia is that one needs a gauge link for the NRQCD matrix element of the 
contribution from the color-octet QQ pair to establish NRQCD factorization beyond one-loop, as shown 
in |20j . We will examine the TMD factorization for p-wave quarkonium in a separate publication. 

To summarize: We have studied the one-loop TMD factorization of 1 5'o-quarkonium production 
in hadron collision at low transverse momentum. We find that the differential cross section can be 
factorized with the TMD gluon distributions, the soft factor and the perturbative coefficient. The TMD 
gluon distributions and the soft factor are consistently defined with QCD operators, the perturbative 
coefficient is determined here at one-loop. In comparison with the factorization derived at tree-level, the 
soft factor is needed to cancel all effects of soft gluons. Our result will be useful not only for extracting 
TMD gluon distributions from experimental data, but also for resumming large log's of gj_ appearing in 
the collinear factorization. 
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